Distinguishing among Scalar Field Models of Dark Energy 



Irit Maor, Ram Brustein 
Department of Physics, Ben-Gurion University, Beer-Sheva 84105, Israel 
e-mai I : irrit , ramybSb gumai I .bgu. ac. il 

We show that various scalar field models of dark energy predict degenerate lu- 

■ minosity distance history of the Universe and thus cannot be distinguished by su- 
O '. 

(•~^ , pernovae measurements alone. In particular, models with a vanishing cosmological 

(N ■ 

constant (the value of the potential at its minimum) are degenerate with models 

• with a positive or negative cosmological constant whose magnitude can be as large 

, as the critical density. Adding information from CMB anisotropy measurements 

^ I does reduce the degeneracy somewhat but not significantly. Our results indicate 

' that a theoretical prior on the preferred form of the potential and the field's initial 

O ■ 

■ conditions may allow to quantitatively estimate model parameters from data. With- 

a^ : 

, out such a theoretical prior only limited qualitative information on the form and 

(N ■ 

I parameters of the potential can be extracted even from very accurate data. 

Ph! I. INTRODUCTION 

^ ■ 



X 



One of the standard methods of interpreting the growing body of evidence from super- 
novae |I| and other measurements that the expansion of the Universe is accelerating, is to 
assume the existence of a dark energy component and to model it using scalar fields (for 
a recent review see 0). This links the expansion history of the Universe to theories of 
fundamental physics. For example, from this perspective the value of the potential at its 
minimum is the cosmological constant (CC). Since at this point there are many theoretical 
ideas about the form of the potential but none that particularly stands out, it would have 
been helpful if the data from cosmological measurements, such as supernovae la, CMB and 
various others, could provide hints about some generic features of the potential. 

Viable scalar field models of dark energy need to have potentials whose energy scale is 
about the critical density ~ 10~^^eV^, and typical scalar field masses about the Hubble mass 
m ~ 10"'^^ eV. In such models typical scalar field variations are about Planck scale rup ~ 10^^ 
GeV, and typical time scales for such variations are about the Hubble time l/i^o ~ 10^^ sec. 
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Whether, and how well, it is possible to determine the parameters and form of potentials 
and the field's dynamical history and future from data beyond such qualitative estimates has 
been addressed previously 0, |, |^, H, 0, ||, ^ |ll| . Weller and Albrecht [Q concluded that 
some potentials could be differentiated using SNAP-like data ||T^. They approximated the 
equation of state (EOS) of each of the models, and showed by likelihood analysis that some 
of the models are distinguishable. Another approach is "reconstruction" |T1|]. Here 

one rewrites the potential as a function of the luminosity-distance {d^) and its derivatives, 
which are in turn functions of the redshift. The potential (and not the EOS) is approximated 
by a fitting function, and statistically tested against a set of accurate dL measurements. The 
efficiency of reconstruction depends on the accuracy of knowing the values of Qm (this is 
needed also when one fits the EOS), and Hq (this is needed only for reconstruction). 

On general grounds we expect that scalar field potentials are less distinguishable than 
their corresponding EOS, because different potentials, with properly adjusted initial condi- 
tions for the field can produce very similar EOS. 

In previous papers [1^, Q we have found that supernovae (SN) measurements are limited 
as a probe of the dark energy EOS wq, due to degeneracies. Specifically, it was shown that 
d^s corresponding to two different wq's are degenerate if both EOS coincide at some point 
at a relatively low red shift, z* (see also []T3| and [T^). The purpose of the present analysis is 
to explore the implications of this degeneracy on the possibility to determine the scalar field 
potential. For a given functional form of potential, we would like to quantify the amount 
by which the parameters of the potential can be varied, and still be indistinguishable by 
accurate SN measurements. Our criteria for indistinguishability between two models is that 
their resulting d^s differ at most by 1% up to redshift z = 2, in accordance with the 
anticipated accuracy of future SN measurements. In addition, we would like to determine 
whether the functional form of the potential can be distinguished or constrained by data. 

We look for degeneracies among potentials using the following procedure: For a given 
class of potentials, we change the parameters as well as the initial values of the field, with 
the constraint that wq at z* remains unchanged. This results in models whose wq cross at 
z*. We know from our previous analysis that in this case the models tend to be degenerate. 
Then we evaluate numerically the differences in the d^s of the models to verify this. 

There are additional sources of degeneracy that we do not consider here. In our procedure 
the value of the potential energy and the value of the kinetic energy remain unchanged. 
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Allowing changes in the potential that are compensated by changes in the initial conditions 
for the kinetic energy will give another dimension of degeneracy. Variation in the value 
of Qjn is yet another degree of freedom, as is relaxing the assumption of a fiat Universe 



and considering the effects of a clumpy Universe Additionally, the value of z* in 

different models can be shifted. We have found that z* varies slowly with the red shift 
depth of the data set, z^ax- The dependence is approximately linear z* = aZmax + P (see 



also |T^), a and (3 are model dependent but a is typically small, about 0.2. And, finally, 



we have loolced only at a class of simple potentials that have two independent parameters. 
Additional parameters in the potential yet again open up new degrees of freedom, each of 
which increases the degeneracy of each of the parameters. Since we have found that this class 
of simple potentials suffers from large degeneracies, we see no phenomenological justification 
for using more complex potentials. If a theoretical prior about the form of the potential and 
initial conditions can be motivated then some of this degeneracy can be removed. 

II. DEGENERACIES OF SCALAR FIELD POTENTIALS 

We consider a fiat Universe with non-relativistic matter (dark matter included) whose 
EOS is Wm = 0, and a dark energy component which we model by a canonically normalized 
and minimally coupled scalar field. Einstein's equations for such a Universe are the following, 

= (p^x' + ix^ifV" + V{(j))) (1) 
xHH' = ^ {p^x^ + x^H^(f)'^) (2) 
x^H^(f)" + {x^HH' - 2xH^)<i)' + ^ = 0, (3) 

where x = 1 + z, = -^2, Pm is the matter energy density today, primes denote derivatives 
with respect to x, and V is the potential of the scalar field. The scalar field's EOS (wq) is 
given by, 

Pq 

Wq = — 
PQ 

_ x^H^cj)'^ - 2V 

~ xm^4>'^ + 2V' ^ ^ 

We choose a model, and vary the parameters of the potential Pi, as well as the values of 
the scalar field at z* keeping its derivative constant, 

cj^iz*) ^ <Piz*) + 5<P (5) 
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V ^ V + 5V{PiJP,, 50) (6) 
[z*) = 6H{z*) = 6pm = 0. (7) 



These variations result in variations in eqs.(^ and (|T]): 

a;2^2^/2 _2{V + 6V) 



Wq + 5wq 



^2^20/2 + 2{y + 5V) 



ff2 = (p„i' + ii'if'/' + (V(4;) + . (9) 

In eq.(II), we have kept terms up to first order in 5V , assuming it is small enough. 

Since we know that if the values of wq{z*) are equal for two models then their luminosity- 
distance history is approximately degenerate, we explore part of the degeneracy in parameter 
space by requiring that 5wq{z*) vanishes. In eq. (H) we ignore higher orders in 5V , so we 
simply require that 5V vanishes to first order. On the other hand, even a small deviation 
from a spatially fiat cosmology will be amplified by the evolution of the solution. We 
therefore need to require that eq. (|) holds exactly. 

The two constraints we are imposing are then 

5V{P,,5P,, (50)« = O (10) 
5V{Pi,5Pi, 50)(^^'^^*) = 0. (11) 

This set of variations and constraints is algebraic and can often be solved analytically. The 
analytical solution connects any given model to a family of associated models, all of which 
have the same wq at z* . The next step is to check numerically how large are the allowed 



parameter variations such that different models in the family are indistinguishable. p4 

Initial conditions are given at z* , but the evolution of each of the models towards z = is 
different, and they end up with a different value of Hq (Hubble parameter at a; = 1, 2; = 0). 
To ensure that all models have the same value of Hq we rescale them. 

Equation (|I|) can be reexpressed in terms of the rescaled variables 

= Qmx^ + \x''H^'^ + V. (13) 
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This means that although initially only the scalar field potential energy is varied, eventually 
all quantities (except Hq) may vary between models. As we shall show later (section [TTTD , it 
turns out that differences in Qm are negligible. This result is not surprising: as was shown in 
(see in particular Fig. 2), this type of degeneracy characterizes models with fixed values 
of flm- Families of models which are degenerate with respect to SN measurements but have 
different values of Qm do not exhibit the enhanced resolution of wq{z*). This means that 
on top of the degeneracy that we will be exploring here, another dimension of degeneracy 
will open up once the uncertainties in Q^n are taken into account. The method described 
here therefore exposes only part of the degeneracy among potentials. 



III. DEGENERACIES AMONG SIMPLE POTENTIALS 



We have analyzed various forms of potentials with two parameters. For all such poten- 
tials there are three independent variations, two parameters and the initial condition for 
0. Enforcing the two constraints eqs. (p!OD,(pAD results in a solution which is a curve in 
a three-dimensional space. Obviously the degeneracy is larger when more parameters are 
allowed. 

One class of potentials we have looked at is V{(j)) = |m^0^ + vq. This is one of the 
standard simple forms of potentials, m is the scalar field mass, Vq is the value of the potential 
at the minimum (the CC). As explained previously we already have an order of magnitude 
estimates for vq and m, but we would like to know whether they can be determined in a 
more quantitative and conclusive way by the data. In particular we would like to know 
whether it is possible to distinguish models with a vanishing CC {vq = 0) from models with 
a non- vanishing CC {vq ^ 0). For this class of quadratic potentials the variation is given 
by: 

1 



V + 6V 



(m + 5m)^(0 + 6(f)y + {vo + 6vo) 



+ 



+ 



1 



-{m6(f) + (pdmY + 5(l)5m{m + 6m){(f) + 6(f)) Sm^l 

2 2 



(14) 



The first bracket is V itself. The second bracket is 5V^^^ which should vanish, and the 
last bracket is SV^^^"''^*^ — SV^^^ which should vanish as well. The solution is a curve in the 
{Svq, 5m, 5(f)) space. On this curve changes in the curvature and the height of the minimum 
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of the potential are compensated by a change in the value of such that the potential energy 
of the field is locally unchanged. 

Figure |I] shows a variety of quadratic potentials on the left, their resulting wq in the 
middle, and the relative difference m on the right. The parameters of the potentials are 
listed in table I. Notice that although the fiducial model has a vanishing CC, vq = 0, it is 
degenerate with models that have Vq of order unity in units of the critical density {rn^Hl). 
The uncertainty in m is of order unity in units of the present value of the Hubble parameter 

All the models have values of Vim which are within ±0.01 from the value of the fiducial. 
This means that we are indeed exploring here the degeneracy due to the integral relation 
between di and wq, and not the degeneracy due to the uncertainty in VLm- As explained 
previously, this is a direct consequence of our method. 

Following the same procedure we have analyzed other popular classes of potentials with 
two parameters. Figure ^ shows similar results for an exponential potential, V = Ae'^"^ . 
For this specific potential only some of the solutions of eqs. (p!OD and (|lT]) can be found 
analytically, hence the degeneracy shown is a smaller than the full degeneracy, as can be 
seen by the reduced range of allowed values of A in table II. Figure ^ shows the results for 
an inverted quadratic potential, V = —^m'^cj)'^ + vq. 

As can be seen from the figures, the deviation in di tends to peak at low redshift, about 
z = 0.4. This results from the following reason: consider two models whose EOS cross 
at z = z*. For redshifts z < z*, the different wq of the models imply different rates of 
expansion. The model whose wq is more negative will expand faster, making the dis of the 
two models diverge. In the range z > z* the trends reverse, making their d^s converge. If 
this were the only source of degeneracy among models, it would have been useful to focus 
measurements in this range of redshifts, z ~ 0.4. Unfortunately a second degeneracy (due 
to the uncertainty in Q^) degrades the extra sensitivity in this region of redshifts. 

So far we have analyzed the uncertainty in determination of parameters under the as- 
sumption that the form of the potential is known. Without a theoretical prior on the form 
of the "true" potential it is important to find whether different classes of potentials can be 
distinguished by the data alone. We will present here only a representative example and not 
attempt a systematic analysis to expose the full degeneracy. 

In Figure ^ we show an example of four different degenerate potentials. All models have 
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FIG. 1: Various quadratic potentials (shown on the left), their resulting wq (middle), and the 
relative difference in (right). The parameters of the potentials are listed in the table I. 





Solid (fiducial) 


Dashed 


Dotted 


Dash-dotted 




1.51 


1.12 


0.98 


1.99 


Vo 
P 





-0.86 


0.20 


0.47 




0.30 


0.31 


0.29 


0.29 



TABLE I: Parameters of quadratic potential plotted in Figure |T|. 
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FIG. 2: Various exponential potentials (shown on the left), their resulting wq (middle), and the 
relative difference in (right). The parameters of the potentials are listed in table II. 
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TABLE II: Parameters of exponential potential plotted in Figure |2|. 
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FIG. 3: Various inverted quadratic potentials (shown on the left), their resulting wq (middle), and 
the relative difference in d-L (right). The parameters of the potentials are listed in table III. 





Solid (fiducial) 


Dashed 


Dotted 


Dash-dotted 




1.51 


1.88 


0.98 


1.08 


Vo 


1.51 


1.09 


1.28 


2.55 




0.30 


0.29 


0.29 


0.31 



TABLE III: Parameters of inverted quadratic potentials plotted in Figure ^. 

Qm = 0.30. The potentials we have used in this example are V = ^m'^cj)'^ + vq (solid), 
V = Ae^^'^ (dotted), V = -\m'^<i? + vq (dashed), and 1/ = a + 60 (dash-dotted). The 
right panel shows that the relative differences in their d^s are less than one percent up to 
z = 2. The left panel shows the potentials. Although the functional forms of the potentials 
are very different and they do not have the same asymptotic behavior, they are nevertheless 
indistinguishable. This can be understood by observing that the patch of potential that was 
probed by the field during the relevant redshift range, (marked in the figure by the shaded 
area) is rather small. In this patch, differences among the various potentials are marginal. 

The fact that a small patch of the potential is actually being probed justifies the use 
of a simple potential with a small number of parameters. A fast rolling field would have 
served as a better probe of the potential since it would have covered a larger patch in a 
given redshift range, but this would have meant a more positive wq which goes against the 
evidence of accelerated expansion. Another possible way to enlarge the size of the probed 
potential patch would have been to measure di over a larger redshift range. The difficulty 
here is that at deeper redshifts the contribution of dark energy to the total energy of the 
Universe is smaller, and its effects become harder to detect. 



9 




FIG. 4: Various potentials (left panel), and their relative differences (right panel). The potentials 
we have used in this example are ^m^c/)^ +vo (solid), Ae~^'^ (dotted), —\m?(tP' +vq (dashed), and 
a + h(j) (dash-dotted). All models have = 0.30. The shaded area in the figure on the left 
indicates the region in which the field moves in the relevant redshift range, < z < 2. 

Our results may seem to disagree with those of where it was shown that different 
forms of potentials could be distinguished using SNAP-like data alone, but we believe that 
they are in fact in full agreement. Our interpretation of the results of is that they 
depend on the choice of specific parameters and initial conditions for each of the potentials 
that makes them distinguishable. This does not mean that for a given set of data, it will 
be possible to determine that a specific class of potentials is preferred over another. It is 
quite easy to construct counter examples by choosing the potentials to look similar in the 
relevant redshift patch. A representative example is shown in Fig. |^. The potentials that 
are used here are the pseudo-Nambu-Goldstone boson potential, V = M'^(cos(0//) + 1) 
(with M'^/irripHl) = 0.82 and f/rrip = 0.6), and the pure exponential, V = Ae~^'^ (with 
A/{mlHl) = 4.43 and rUpB = 1.78). 

IV. CONSIDERING CMB 

An additional measurement that can possibly probe the dark energy EOS and the scalar 
field potential is that of the cosmic microwave background (CMB) anisotropy. This mea- 
surement is expected to be improved with the upcoming MAP and PLANCK missions. 

Using the CMB measurements to probe scalar field models requires additional theoretical 
assumptions about the evolution in the range from z ~ 2 to z ~ 1000. For example, that 
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FIG. 5: Two potentials (shown on the left), and their relative (11 differences (right). The potentials 
we have used in this example are v = M^(cos(|) + 1) (solid), and Ae'^'l' (dashed). Both models 
have 0,m = 0.30. The shaded area in the left panel indicates the region in which the field moves in 
the relevant redshift range, < z < 2. 



the potentials are "tracker" potentials |T9|, or that the scalar field has negligible influence 
above some redshift (say, z = 2). 

It is well-known that the EOS cannot be determined using CMB anisotropy data alone due 



to degeneracies. A thorough analysis was done in |2^, who concluded that this degeneracy 
persists after combining CMB and SN measurements. In [|14| it was shown that the inclusion 
of CMB constraints in the analysis of the dark energy EOS does improve the resolution 
somewhat, but not significantly. 

In a simple minded estimate here, we would like to show that under reasonable assump- 
tions, CMB measurements will not help to significantly constrain the potential either. Rather 
than performing an extensive numerical search, as in [| 



Plj |2^ , we use a different strategy: 
we treat CMB measurements as effectively providing one additional point on the Hubble 
diagram at the last scattering redshift zis ~ 1000. The method cannot provide numerically 
accurate results, but it does highlight the theoretical degeneracy that limits the ability of 
CMB measurements to constrain the parameters and functional form of the potential. 

Treating the CMB measurements as effectively one additional point on the Hubble dia- 
gram at redshift z ~ 1000 can be implemented in a simple way because the angular distance 
and the luminosity distance are related to each other = (1 + z^dA- The observed angular 
size of any feature in the CMB, 6, is related to its physical size d by the angular distance 
c?A to last scattering surface, 9dA = d. The CMB spectrum yields a series of acoustic peaks, 
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located at /„ = ^dA, n being an integer, and 5* is the sound horizon at last scattering 
surface. We are going to ignore uncertainties in the sound horizon S, which depends on the 
composition of the universe at last scattering surface. 

To estimate the measurement accuracy of the CMB point on the Hubble diagram we 
may treat the position of the first peak /i, as a direct measurement of the angular distance, 
therefore the relative errors in li and in the angular distance are equal ^ = The 
position of the first peak is already measured rather well, to about 3% at the la level (see. 



for example, |23|). This is expected to be somewhat improved by MAP and PLANCK, so 
we may expect eventually a 3a error in the few percent range. 



Another way of estimating the accuracy of the CMB point is to use the results of [14 



where it was shown that for models which have a linear wq for z < 2 and wq = const, from 
z = 2 to last scattering surface, a derivative as large as t^i = ±| could not be distinguished 
at the 3a level, assuming full sky coverage and cosmic variance limited measurement. This 
means that the most accurate CMB measurements do not distinguish models which lead to 
a difference of (+1.6%, -3.2%) in their dis. Thus a 3a measurement error estimate of a few 
percent in di seems reasonable. 

According to the preceding discussion we may define models to be indistinguishable by the 
CMB if their di difference is less than a few percent at zig. We may check now whether this 
imposes additional constraints on our models. We find that it does, but not to a significant 
degree. 

Since our toy models are not strictly chosen as tracker models, evolving the equations of 
motion backwards from 2; ~ 0.3 up to 2; ~ 1000 gives nonsensical results, the dark energy 
typically becomes dominant, with wq approaching 1. In order to to avoid such undesirable 
and observationally excluded behavior, the potential needs to be changed in the regions that 
are probed by the field in the redshift range z ^ 1 such that the kinetic energy remains 
acceptably small, so that wq remains negative and Qq becomes sub-dominant. 

Instead of building the potential piece-wise, we have chosen another approach in the spirit 
of [|l^: we let the equations of motion evolve until z = 2 (the redshift range relevant to SN 
measurement), and for z > 2 we put by hand wq = for all models. We then calculate 
numerically the relative differences in d^ at z = 1000. Our cutoff procedure changes the 
z-dependence of wq at a value of z where the dark energy density is negligible, which makes 
the details of the ^-dependent cutoff unimportant. We expect other cutoff procedures to 
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give similar results. 

We find that the relative differences in di for our models are within a few percent, in 
agreement with the semi-analytic argument that we present shortly. For the potentials 
shown in Fig. ^, we find = —0.9% for the exponential potential, —3.4% for the inverted 
quadratic potential, and —1.6% for the linear potential. As explained previously this is just 
at the limit of CMB resolution. For the quadratic potential examples in Fig. |I] we find 
the following di differences from the fiducial, dashed: 2.5%, dot-dashed: —6.2%, dotted: 
8.0%, for the exponential potentials of Fig. || we find: dashed: 4.5%, dotted: 6.8%, and 
for the inverted quadratic potential of Fig. ^ we find dashed: —7.8%, dotted: —8.1%, and 
dot-dashed: 5.1%. 

For most of the models we see that CMB is expected to further constrain the parameters if 
the form of the potential is known. The larger d^ differences are for models whose differences 
at 2; ~ 2 are larger, in agreement with the following semi- analytic argument. Given the 
approximate nature of our measurement accuracy estimate, the large number of additional 
sources of degeneracy that we have neglected, and the additional theoretical assumptions 
that go into the analysis with the CMB point added it is not possible, and we believe that it is 
not necessary to determine in a more quantitative way the amount by which the degeneracy 
is reduced. 

We may estimate in a rough semi-analytic way the relative errors in d^ at z ^ 1000 as 
follows. Luminosity distance at 1 + zis = xis ~ 1000 for models with wq{x > 3) = is given 
by El, 

I f^if dx ( dy \ 

dLixis) = xis^jl + g U + ea;p(3 

= xis.JlT~gJ^ (^g + exp{3 w;q^)^ 
/ /"^'s dx ( dy \ 

= ^dL{3)+xis^lT~gJ^"^(^g + exp{3j^ w;q^)^ 

1 /2 

= ^dL{3) + 2xu I ^ I , (15) 

3 ^' \V3 yWs) [g + exp{3j^WQ'f)) ^ ^ 

where g is the current ratio of matter to dark energy densities, and we have set Hq to unity. 
The second term can be expressed in terms of the luminosity distance at a; = 3 and its 



13 



derivative: 

1+g \ ^3/2 d di 



+ exp{3 WQ^) J dx 



^'^'i: - ^ (16) 



X 



x=3 



dL{xis) ^ - + 2xi,rf'j3). (17) 

Now we would like to examine the difference in diixis) for two models. For a simple 
and rough error estimate, we may use c?^(3) = Ci^^^ and A(i'j^(3) = C2^^y^ with ci, C2 
numerical coefficients of order unity. So finally we obtain 



/ AdA 2c2 - 1 / AdL 



(18) 



\ di J 2ci — 1 \ dL I „ 

Note that the relative di error at xis is independent of xis- Since we consider models whose 
maximal relative error at z = 2 is one percent, they will have a relative error in di of about 
a few percent at zis, depending on the values of ci and C2. This estimate agrees with the 
numerical examples that we have listed above. 

We conclude that once physically reasonable constraints are imposed on the potential, 
CMB measurements do not significantly constrain the parameter space of the potentials, 
although they are expected to reduce it somewhat. Our simple minded analysis strengthens 
the case first made in [EO] where it was shown that only some average EOS can be measured. 



and agrees with [118], E2 



V. CONCLUSIONS 



We have found that it is not possible to obtain precise quantitative estimates for pa- 
rameters of scalar field models of dark energy from data alone beyond the obvious order of 
magnitude estimates. This is due to theoretical degeneracies, which would persist even with 
expected future data from the most accurate SN and CMB measurements. 

Theoretical prior knowledge or assumptions on the form of the potential and the field's 
initial conditions (preferably leaving a total of just two free parameters) may allow a more 
quantitative determination. For example, assuming that the field is at rest at the bottom of 
the potential is equivalent to having a pure cosmological constant. In this case the magnitude 
of the cosmological constant can be determined with accurate data to within a few percent. 
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